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Abstract 

A new notion of stochastic ordering is introduced to compare multi- 
variate stochastic risk models with respect to extreme portfolio losses. 
In the framework of multivariate regular variation comparison criteria 
are derived in terms of ordering conditions on the spectral measures, 
which allows for analytical or numerical verification in practical appli- 
cations. Additional comparison criteria in terms of further stochastic 
orderings are derived. The application examples include worst case 
and best case scenarios, elliptically contoured distributions, and mul- 
tivariate regularly varying models with Gumbel, Archimedean, and 
Galambos copulas. 

1 Introduction 

This paper is dedicated to the comparison of multivariate probability distri- 
butions with respect to extreme portfolio losses. A new notion of stochastic 
ordering named asymptotic portfolio loss order (^ ap i) is introduced. Specially 
designed for the ordering of stochastic risk models with respect to extreme 
portfolio losses, this notion allows to compare the inherent extreme portfolio 
risks associated with different model parameters such as correlations, other 
kinds of dependence coefficients, or diffusion parameters. 

In a recent paper of Mainik and Riischendorf (2010) the notion of extreme 
risk index has been introduced in the framework of multivariate regular vari- 
ation. This index, denoted by 7^, is a functional of the vector £ of portfolio 
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weights and of the characteristics of the multivariate regular variation of 
X given by the tail index a and the spectral measure VP. It measures the 
sensitivity of the portfolio loss to extremal events and characterizes the prob- 
ability distribution of extreme losses. In particular, it serves to determine the 
optimal portfolio diversification with respect to extreme losses. Within the 
framework of multivariate regular variation the notion of asymptotic portfolio 
loss ordering introduced in this paper is tightly related to model comparison 
in terms of the extreme risk index 7^. Thus this paper can be seen as a 
supplement of the previous one, allowing to order multivariate risk models 
with respect to their extremal portfolio loss behaviour. 

In Section 2 of the present paper we introduce the asymptotic portfolio 
loss order ^ ap i and highlight some relationships to further well-known order- 
ing notions. It turns out that even strong dependence and convexity orders do 
not imply the asymptotic portfolio loss order in general. We present counter- 
examples, based on the the inversion of diversification effects in models with 
infinite loss expectations. Another example of particular interest discussed 
here is given by the elliptical distributions. In this model family we estab- 
lish a precise criterion for the asymptotic portfolio loss order, which perfectly 
accords with the classical results upon other well-known order relations. Sec- 
tion 3 is devoted to multivariate regularly varying models. We discuss the 
relationship between the asymptotic portfolio loss order and the comparison 
of the extreme risk index and characterize ^ ap i in terms of a suitable or- 
dering of the canonical spectral measures. These findings allow to establish 
sufficient conditions for ;< ap i in terms of spectral measures, which can be ver- 
ified by analytical or numerical methods. In particular, we characterize the 
dependence structures that yield the best and the worst possible diversifica- 
tion effects for a multivariate regularly varying risk vector X in with tail 
index a. For a > 1 the best case is given by the asymptotic independence 
and the worst case is the asymptotic comonotonicity. The result for a < 1 is 
exactly the opposite (cf. Theorem 3.7 and Corollary 3.8). Restricting X to 

means that X represents only the losses, whereas the gains are modelled 
separately. This modelling approach is particularly suitable for applications 
in insurance, operational risk, and credit risk. If X represents both losses 
and gains, these results remain valid if the extremal behaviour of the gains 
is weaker than that of the losses, so that there is no loss-gain compensation 
for extremal events. In Section 4 we discuss the interconnections between 
^ ap i or ordered canonical spectral measures and other well-known notions 
of stochastic ordering. Ordering of canonical spectral measures allows to 
conclude ^ ap i from the (directionally) convex or the supermodular order. It 
is not obvious how to obtain this implication in a general setting. Finally, 
in Section 5 we present a series of examples with graphics illustrating the 
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numerical results upon the ordering of spectral measures. The relationship 
to spectral measures provides a useful numerical tool to establish ^ apl in 
practical applications. 



2 Asymptotic portfolio loss ordering 

To compare stochastic risk models with respect to extreme portfolio losses, 
we introduce the asymptotic portfolio loss order z< ap i. This order relation 
is designed for the analysis of the asymptotic diversification effects and the 
identification of models that generate portfolio risks with stronger extremal 
behaviour. 

Before stating the definition, some basic notation is needed. Focusing 
on risks, let X be a random loss vector with values in M. d , i.e., let positive 
values of the components i = \,...,d, represent losses and let negative 
values of represent gains of some risky assets. Following the intuition of 
diversifying a unit capital over several assets, we restrict the set of portfolios 
to the unit simplex in IR d : 

The portfolio loss resulting from a random vector X and the portfolio £ is 
given by the scalar product of £ and X. In the sequel it will be denoted by 

ex. 

Definition 2.1. Let X and Y be d- dimensional random vectors. Then X is 
called smaller than Y in asymptotic portfolio loss order, X ^ ap i Y , if 

V£eE d limsup ^ f ! T 5 > ? <!■ (1) 
Here, jj is defined to be 1. 

Remark 2.2. (a) Although designed for random vectors, ^ ap i is also defined 
for random variables. In this case, the portfolio set has only one element, 
£! = {!}. 

(b) It is obvious that ^ ap i is invariant under componentwise rescaling. Let 
vx denote the componentwise product of v,x e M. d : 

(AW,...,/'^), (2) 

Then it is easy to see that X ^ ap i Y implies vX ^ ap i vY for all v G R+. 
Hence condition (1) can be equivalently stated for £ e M+. 
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The ordering statement X ^ apl Y means that for all portfolios the 
portfolio loss £ T X is asymptotically smaller C, T Y. Thus ;< ap i concerns only 
the extreme portfolio losses. In consequence, this order relation is weaker 
than the (usual) stochastic ordering ^ st of the portfolio losses: 

£ T X =<* t T Y for all CeE d implies X ^ apl Y. (3) 

Here, for real random variables U, V the stochastic ordering U ^ st V is 
defined by 

Vt G R P{U >t}< P{V > t}. (4) 

Some related, well-known stochastic orderings (cf. Miiller and Stoyan, 
2002; Shaked and Shanthikumar, 1997) are collected in the following list. 
Remind that / : R d — > R is called supermodular if 

Vx,yeR d f(xAy) + f(xVy)>f(x) + f(y). (5) 

Definition 2.3. Let X , Y be random vectors in R d . Then X is said to be 
smaller than Y in 

(a) (increasing) convex order, X ^ cx Y (X ^ icx Y), ifEf(X) < Ef(Y) for 
all (increasing) convex functions f : R d h-> R such that the expectations 
exist; 

(b) linear convex order, X ^ kx Y , if £ T X ^ cx £ T Y for all £ G R°7 

(cj positive linear convex order, X ^ p i cx Y, if i T X ^ cx £ T F for all ( G ffj.; 

(d) supermodular order X ;< sm Y, ifEf(X) < Ef(Y) for all supermodular 
functions f : M d — > R sfic/i that the expectations exist; 

(e) directionally convex order, X ^dcx Y , ifEf(X) < Ef(Y) for all direc- 
tionally convex, i.e., supermodular and componentwise convex functions 
f : R d — > R such that the expectations exist. 

The stochastic orderings listed in Definition 2.3 are useful for describing 
the risk induced by larger diffusion (convex risk) as well as the risk induced by 
positive dependence (supermodular and directionally convex). The following 
implications are known to hold generally for random vectors X, Y in R d : 

(a) (X ^ sm Y) (X ^ dcx Y) (X ^ plcx Y) 

(b) (X ^ cx Y) (X ^ lcx Y) =S> (X ^ plcx Y) 

(c) (X ^ icx Y) (X ^ plcx Y) 
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Remark 2.4. (a) It is easy to see that the usual stochastic order ^ st implies 
^ ap i in the univariate case. 

(b) In spite of being strong risk comparison orders, the order relations out- 
lined in Definition 2.3 do not imply ;< ap i in general. For instance, it 
is known that the comonotonic dependence structure is the worst case 
with respect to the strong supermodular ordering ^ sm , whereas it is not 
necessarily the worst case with respect to ^ ap i (cf. Examples 5.1 and 5.2). 

The following proposition helps to establish sufficient criteria for ^ ap i 
in the univariate case. To obtain multivariate results, it can be separately 
applied to each portfolio loss £ T X for £ e S d . 

Proposition 2.5. Let R 2 > be real random variables and let V be a 
real random variable independent of Ri, % = 1,2. 

(a) If R\ ^ ap i R 2 an d V < K for some constant K, then 

RiV ^ apl R 2 V. (6) 

(b) IfRi dst Ri, then 

(RiV) + dst and (R 2 V)_ dst (RiV)- (7) 

In addition, ifV and Ri are integrable and EV > 0, then 

RiV ^ icx R 2 V. (8) 
Moreover, if EV = 0, then R{V ^ cx R 2 V . 
Proof. 

Part (a). Since RiV ^ ap i R 2 V is trivial for V < 0, we assume that 
P{V > 0} > 0. Hence V < K implies for all t > 



P{R 1 V>t}=[ P{i?x > t/v}dP v (v) 

= [ f(t/v)P{R 2 >t/v}dP v (v), (9) 

J(0,K) 



where 

I[z) -p{R 2>z y 

An obvious consequence of (9) is the inequality 

P{RiV >t}< sup {f\z) : z > t/K} ■ P {R 2 V > t} (10) 
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Since Ri ^ ap i R 2 is equivalent to limsup 2 _ >00 f(z) < 1, we obtain 

hmsup — r- < 1. 

t^J P{R 2 V >t}~ 

Part (b). By the well-known coupling principle for the stochastic ordering 
^ st we may assume without loss of generality that Ri < R 2 pointwise on the 
underlying probability space. This implies 

P{i?i\/ > t} < P{R 2 V >t}, t> 0, 

and, similarly, 

P{RiV <t}< P{R 2 V <t}, t< 0. 

In consequence we obtain (7). 

From the proof of (7) it follows that the distribution functions of the 
products RiV, i = 1,2, satisfy the cut criterion of Karlin-Novikov (cf. 
Shaked and Shanthikumar, 1994, Theorem 2. A. 17 and Miiller and Stoyan, 
2002, Theorem 1.5.17) Hence we obtain 

RiV^ icx R 2 V. (11) 

If EV = 0, then E^V] = E[R 2 V] and therefore 

RiV ^ cx R 2 V. (12) 

□ 

Remark 2.6. (a) Note that (7) implies (without assuming the existence of 
moments) that (R 2 V) + ^dccx (RiV) + where ^decx denotes the decreasing 
convex order. Similarly one obtains (R 2 V)- ^i CX (RiV)- 

(b) If f(t) := P{Ri > t}/P{R 2 > t} < C < 00 and R 1 ^ api R 2 , then 
RiV ^ ap i R 2 V. 

(c) A related problem is the ordering of products RVi for R > with V\ 
and V 2 independent of R. In the special case when R is regularly varying 
with tail index a > 0, i.e., 

lim Pj^>te} = 

t^oo P{R >t} ' ' V ; 

exact criteria for ^ ap i can be obtained from Breiman's Theorem (cf. 
Resnick, 2007, Proposition 7.5). If E(\/;)+ +£ < 00 for i = 1,2 and some 
e > 0, then 

t^oo P{ J R>t} LV * ;+J 
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This yields 

g > tj = e [ova 

t^oc P{RV 2 >t} E [(V 2 )%] ' 

An important class of stochastic models with various applications are el- 
liptical distributions, which are natural generalizations of multivariate normal 
distributions. A random vector X G M, d is called elliptically distributed, if 
there exist \i G ~R d and a dx d matrix A such that X has a representation of 
the form 

X = fi + RAU, (14) 
where U is uniformly distributed on the Euclidean unit sphere §2, 

S d 2 = {xeR d : \\x\\ 2 = 1} , 

and R is a non-negative random variable independent of U . By definition we 
have 

E\\X\\ 2 2 < 00 0ER 2 < 00, (15) 

and in this case 

Cov(X) = Vax(R)AA T . (16) 

The matrix C := AA T is unique except for a constant factor and is also called 
the generalized covariance matrix of X. We denote the elliptical distribution 
constructed according to (14) by £(//, C, Fr), where Fr is the distribution of 
R. 

A classical stochastic ordering result going back to Anderson (1955) and 
Fefferman et al. (1972) (cf. Tong, 1980, p. 70) says that positive semidefinite 
ordering of the generalized covariance matrices C\ -< w a C 2 , defined as 

v£ g M. d e^^ec^, (17) 

implies symmetric convex ordering if the location parameter \x and the dis- 
tribution Fr of the radial factor are fixed: 

£{[/,, Ci, Fr) ^sjmma S(ijl,C 2 ,F r ). (18) 

It is also known that for elliptical random vectors X ~ £(/i, C, Fr) the multi- 
variate distribution function F(x) := P{Xl < xi, . . . ,X d < x d } is increasing 
in Cij for i 7^ j, where C = (C it j) (see, e.g., Joe, 1997, Theorem 2.21). 

The following result is concerned with the asymptotic portfolio loss or- 
dering ^ ap i for elliptical distributions. 
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Theorem 2.7. Let X = pi + RxA-JJ , Y = p 2 + R2A2U be elliptically dis- 
tributed with generalized covariances Ci := AiAj . If 

p 1 < p 2 , R\ ^a P i R2, (19) 

and 

v^es d ec 1 z<ec 2 z, (20) 

then 

X ^ apl Y. (21) 

Proof. It suffices to show that £ T Y ;< ap i ^ T Y for an arbitrary portfolio 
£ G S d . Furthermore, without loss of generality we can assume p\ = p2 = 0. 
For i — 1,2 and £ G £ d denote 



a, = a t (0 := (£ T a£) 

and 



1/2 



Then, by definition of elliptical distributions, we have 

CX = R iaiVl U and £ T Y = R 2 a 2 v 2 U. (22) 

Since the vectors Vi = i>i(£) have unit length by construction, the random 
variables vJJ are orthogonal projections of U ~ unif^f) on vectors of unit 
length. Symmetry arguments yield that the distribution of UjC/ is indepen- 
dent of vt and that v,JJ = (1,0,..., 0) T U = U (1 l 
Thus we have 

^X^a^V and £ T Y = a 2 R 2 V 

with V := By assumption we have a\ < 02 and i? x ;< ap i R 2 . Applying 

Proposition 2.5(a) we obtain £ T X ;< apl £ T F. □ 

Remark 2.8. (a) It should be noted that condition (20) is indeed weaker than 
(17). Let — 1 < pi < p 2 < 1 and consider covariance matrices 

(':= I 1 P ! 1. / 1.2 



Pi 1 

Straightforward calculations show that C; satisfy (20), but not (17). 

(b) For subexponentially distributed Ri the assumption p\ < p 2 in (19) can 
be omitted. 
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3 Multivariate regular variation: ^ ap i in terms 
of spectral measures 

This section is concerned with the characterization of the asymptotic port- 
folio loss order ^ ap i in the framework of multivariate regular variation. The 
results obtained here highlight the influence of the tail index a and the spec- 
tral measure \I/ on z< ap i, with primary focus put on dependence structures 
captured by It is shown that ^ ap i corresponds to a family of order re- 
lations on the set of canonical spectral measures and that these order re- 
lations are intimately related to the extreme risk index 7^ introduced in 
Mainik and Ruschendorf (2010) and Mainik (2010). 

The main result of this section is stated in Theorem 3.6, providing cri- 
teria for X ^ ap i Y in terms of componentwise ordering -< apl yO) f or 
i — 1, . . . , d and ordering of canonical spectral measures. A particular conse- 
quence of these criteria is the characterization of the dependence structures 
that yield the best and the worst possible diversification effects for random 
vectors in (cf. Theorem 3.7 and Corollary 3.8). Another application con- 
cerns elliptical distributions. Combining Theorem 3.6 with results on ^ ap i 
obtained in Theorem 2.7, we obtain ordering of the corresponding canonical 
spectral measures. 

Recall the notions of regular variation. In the univariate case it can be 
defined separately for the lower and the upper tail of a random variable 
via (13). A random vector X taking values in R d is called multivariate 
regularly varying with tail index a G (0, 00) if there exist a sequence a n — > 00 
and a (non-zero) Radon measure v on the Borel a-field #([—00, oo] rf \ {0}) 
such that 00, oo] d \ M. d ) = and, as n — > 00, 

nP"" lx Ai/onB([-oo,oo] d \{0}), (23) 

where -4- denotes the vague convergence of Radon measures and P a ™ x is the 
probability distribution of a~ l X. 

It should be noted that random vectors with non-negative components 
yield limit measures v that are concentrated on [0, oo] d \ {0}. Therefore 
multivariate regular variation in this special case can also be defined by 
vague convergence on £>([0, oo] d \ {0}). 

Many popular distribution models are multivariate regularly varying. In 
particular, according to Hult and Lindskog (2002), multivariate regular vari- 
ation of an elliptical distribution S(fi,C, F R ) is equivalent to the regular 
variation of the radial factor R and the tail index a is inherited from R. 
Other popular examples are obtained by endowing regularly varying margins 



9 



X® with an appropriate copula Wiithrich (cf. 2003); Alink et al. (cf. 2004); 
Barbe et al. (cf. 2006) 

For a full account of technical details related to the notion of multivariate 
regular variation, vague convergence, and the Borel a-fields on the punctured 
spaces [— oo, oo] d \ {0} and [0, oo] d \ {0} the reader is referred to Resnick 
(2007). 

It is well known that the limit measure v obtained in (23) is unique 
except for a constant factor, has a singularity in the origin in the sense that 
u((—e,e) d ) = oo for any e > 0, and exhibits the scaling property 

v{tA) = r a v{A) (24) 

for all sets A e B ([— oo, oo] d \ {0}) that are bounded away from 0. 

It is also well known that (23) implies that the random variable ||X|| with 
an arbitrary norm ||-|| on M. d is univariate regularly varying with tail index 
a. Moreover, the sequence a n can always be chosen as 

a n :=Ffc {l (l-l/n), (25) 

where ^jxy i s the quantile function of \\X\\. The resulting limit measure v is 
normalized on the set A\\.\\ := {i G K d : ||x|| > 1} by 

V (j4||.||) = 1. (26) 

Thus, after normalizing v by (26), the scaling relation (24) yields an 
equivalent rewriting of the multivariate regular variation condition (23) in 
terms of weak convergence: 

C {t^X | ||X|| > t) A v\ Am on B (A H ) (27) 

for t — > oo, where ^U| M | is the restriction of v to the set A\\.\\. 

Additionally to (23) it is assumed that the limit measure v is non-degen- 
erate in the following sense: 

v({x ER d : \x^\> l}) > 0, i = l,...,d. (28) 

This assumption ensures that all asset losses are relevant for the extremes 
of the portfolio loss £ T X. If (28) is satisfied in the upper tail region, i.e., if 

v({x G R d : x {i) > 1}) > 0, i = l,...,d, (29) 

then v also characterizes the asymptotic distribution of the componentwise 
maxima M n := {M {1 \ . . . , M«) with M« := m&x{x[ i] , . . . , X^} by the 
limit relation 

P {a- l M n E [~oo,x]} A exp (-v ([-oo, oo\ d \ [-oo, x})) (30) 
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for x G (0, oo] d . Therefore v is called exponent measure. For more details 
concerning the asymptotic distributions of maxima the reader is referred 
to Resnick (1987) and de Haan and Ferreira (2006). 

Another consequence of the scaling property (24) is the product repre- 
sentation of v in polar coordinates 

(r, s) := t(x) := (||a;||, ||a;|| _1 :r) 

with respect to an arbitrary norm ||-|| on R d . The induced measure v T := 
v o t -1 necessarily satisfies 

v T = c-p a ®^ (31) 

with the constant factor 

c = v(A H ) > 0, 
the measure p a on (0, oo] defined by 

p a {{x,oo]) := x~ a , x£(0,oo], (32) 

and a probability measure \1/ on the unit sphere Sjf ,, with respect to ||-||, 

Sf.|| := {sGK d : ||s|| = 1} . 

The measure ^ is called spectral measure of v or X . Since the term "spectral 
measure" is already used in other areas, \1/ is also referred to as angular 
measure. In the special case of M^-valued random vectors X it may be 
convenient to reduce the domain of ^ to fl IR^. 

Although the domain of the spectral measure \I/ depends on the norm ||-|| 
underlying the polar coordinates, the representation (31) is norm-independent 
in the following sense: if (31) holds for some norm ||-||, then it also holds for 
any other norm ||-|| that is equivalent to ||-||. The tail index a is the same 
and the spectral measure \l/ on the unit sphere corresponding to ||-|| is 
obtained from \1/ by the following transformation: 

V = V T , T(s) := WsW^s. 

Finally, it should be noted that multivariate regular variation of the loss 
vector X is intimately related with the univariate regular variation of port- 
folio losses £ T X. As shown in Basrak et al. (2002), multivariate regular vari- 
ation of X implies existence of a portfolio vector £ £ ^ d such that £jX is 
regularly varying with tail index a and any portfolio loss £ T X satisfies 

lim P U X >t | = u £ ) e [o oo). (33) 
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This means that all portfolio losses £ T X are either regularly varying with 
tail index a or asymptotically negligible compared to £jX. 

Moreover, it is also worth a remark that for Revalued random vectors X 
the converse implication is true in the sense that (33) and univariate regular 
variation of ( I imply multivariate regular variation of the random vector 
X. This sort of Cramer- Wold theorem was established in Basrak et al. (2002) 
and Boman and Lindskog (2009). 

Under the assumption of multivariate regular variation of X the extreme 
risk index 7^ = J^(X) is defined as 

-*w=felSK£ir (34) 

In Mainik and Riischendorf (2010) the random vector X is restricted to M.+ 
and the portfolio vector £ is restricted to S d . The general case with X in M. d 
and possible negative portfolio weights, i.e., short positions, is considered in 
Mainik (2010). Normalizing the exponent measure v by (26), one obtains 

#)=.({iGK d :( T i>l}). (35) 

Rewriting this representation in terms of the spectral measure \P and the tail 
index a yields 

7€= / (Cs) a + dns). (36) 

j si 

Denoting the integrand by f^ a , we will write this representation as 7^ = 

The extreme risk index 7^ (X) allows to compare the risk of different 
portfolios. It is easy to see that (34) implies 

P&X > t} _ jOXl (3 7) 

t l ™p{£x>t}-^ 2 (xy (37) 

Thus, by construction, ordering of the extreme risk index 7^ is related to the 
asymptotic portfolio loss order z< ap i- 

However, designed for the comparison of different portfolio risks within 
one model, the extreme risk index 7^ cannot be directly applied to the com- 
parison of different models. The major problem is the standardization by 
P{||X||i > t} in (34). Indeed, since P{||X||i > t} also depends on the spec- 
tral measure of X, criteria for ^ ap i in terms of 7^ demand the specification 
of the limit 

lim^fcO. 

*->«, p{||y|| 1 > t} 
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Another technical issue arises from the invariance of ^ ap i under com- 
ponentwise rescalings. Since the spectral measure \& does not exhibit this 
property, ordering of spectral measures needs additional normalization of 
margins that makes it consistent with z^ ap i. To solve these problems, we use 
an alternative representation of 7^ in terms of the so-called canonical spectral 
measure \f r *, which has standardized marginal weights. 

This representation is closely related to the asymptotic risk aggregation 
coefficient discussed by Barbe et al. (2006). Furthermore, the link between 
the canonical spectral measure and extreme value copulas allows to transfer 
ordering results for copulas into the ^ ap i setting. These results are presented 
in Section 4. 

To reduce the problem to the essentials, we start with the observation 
that ^ ap i is trivial for multivariate regularly varying random vectors with 
different tail indices and non-degenerate portfolio losses. 

Proposition 3.1. Let X and Y be multivariate regularly varying on M. d and 
assume that 7^ (Y) > for all £ e S d . 



(a) V 



hg^U, (38) 

t^oo p{||y||i > *} ' v ; 



then X ;< apl Y. 

(b) If ax > oly, then X ^ apl Y . 
Proof. 

(a) Using relation (34) we obtain 

p {Cx > t} 

h T™ p ?{CY>t } 

= hmsup ^{^>^ . P{imii>t} . P{||X|| 1 >t} 

^^Pdixn^o P{CY>t} p{||y||i>*} 

75 on ^oo p p{iiriii>t} 
= 0. 

(b) Recall that multivariate regular variation of X implies regular variation 
of || X ||i with tail index ax- Analogously, ||V||i is regularly varying with 
tail index ay- Finally, ax > ay yields (38) and by (3.1) we obtain 

A^ apl r. □ 
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Thus the primary setting for studying the influence of dependence struc- 
tures on the ordering of extreme portfolio losses is the case of random vari- 
ables X and Y with equal tail indices: 

ax = CHy —'■ 01. 

In the framework of multivariate regular variation, asymptotic dependence 
in the tail region is characterized by the spectral measure ^ or its canonical 
version The canonical exponent measure v* of X is obtained from the 
exponent measure v as 

v* = v o T 

with the transformation T : M. d — > M. d defined by 

T(x) := (T a • , • • • , T a (u(B d ) ■ , (39) 

where 

T a (t) := (t]l a - and B, := {x E R d : \x®\ > 1} . (40) 

Furthermore, v* exhibits the scaling property 

u*(tA) = r 1 u*(A), t>0, 
and, analogously to (31), has a product structure in polar coordinates: 

v* o r _1 = pi <g> (41) 

The measure is the canonical spectral measure of X. 

Since ^ ap i and ^* are invariant under componentwise rescalings, the 
canonical spectral measure VP* is more suitable for the characterization of 
i^api- The following lemma provides a representation of the extreme risk 
index 7^ in terms of Note that the formulation makes use of the compo- 
nentwise product notation (2). 

Proposition 3.2. Let X be multivariate regularly varying on M, d with tail 
index a G (0, 00). If X satisfies the non- degeneracy condition (28), then 

7 € pQ= / gs, a (vs)d**(s), (42) 

where ty* denotes the canonical spectral measure of X , the rescaling vector 
v = (v^\ . . . ,v^) is defined by 

v^:=( let (X) + 1 ^(X)), (43) 

and the function g^ a : R d — > R is defined as 

^):=(E^-((^)r-(^- a )) • ( 44 ) 
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Proof. Denote := {i G l d : ( T i > 1}. Then, by definition of v\ 

7c (X) = u(A^) 

= v* {x G R d : T(ar) G A u } 

= // 1 (e T T(rs) > 1} d Pl (r)d^*(s). (45) 



It is easy to see that (40) implies T a (rt) = r 1/a T a (t) for r > and t G R. 
Consequently, (39) yields 

T(rx) = r 1/a T(x) (46) 
for r > and x G lR d . Applying (46) to (45), one obtains 

7«P0 = / / 1 (r 1/a e T T( S ) > 1} d Pl (r) d**( S ) 

</Sf </(0,oo) 

= / / i {ens) > o} i {r > (ens)y a \ d Pl (r)d**( S ) 

JS( J(0,oo) 1 J 

= [ i{ens)>o} {fT( S )) a d**( S ) 

Jsf 

= [ (eT(s))l dV(s). (47) 
Finally, consider the sets Bi defined in (40). It is easy to see that 

Hence 

= !k..: (VS) . □ 



As already mentioned above, ^ ap i and are invariant under rescaling 
of components. Consequently, characterization of ^ ap i can be reduced to 
the case when the marginal weights = 7 e .(X) + r y_ e .(X) in (42) are 
standardized by 

!' S£ p { { H>|} = l (48) 
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This condition will be referred to as the balanced tails condition. The fol- 
lowing result shows that this condition significantly simplifies the represen- 
tation (42). 

Proposition 3.3. Suppose that X is multivariate regularly varying on M. d 
with tail index a G (0, oo). 

(a) If X has balanced tails in the sense of (48), then 

7*P0 



7ei pO + 7 _ ei pO 



(49) 



(b) The non- degeneracy condition (28) is equivalent to the existence of a 
vector w G (0, oo) d such that wX has balanced tails. 

(c) The extreme risk index 75 of the rescaled vector wX obtained in part (b) 
satisfies 

Wrf)+7- t ,(rf) = ***-■ (50) 
Proof. Part (a). Consider the integrand g^ a (vs) in the representation (42): 

9*A VS ) = ■ ((u (i) s (i) ) + a - 

The balanced tails condition (48) implies that X is non-degenerate in the 
sense of (28). Furthermore, all weights in the representation (42) are 
equal: 

1 , P{\X^\>t}/P{\\X\\ 1 >t} 7ei (X)+ 7 _ e ,(X) 

tiJS, P {\XU)\ >t}/P{\\X\\ 1 >t} 7ej (X) + 7 _ ej (X) 
v {i) 

Hence g^ a (vs) simplifies to 

9tA vs ) = u (1) #,a(s) 

= ( 7ei (X)+ 7 „ ei (X))^ a ( S ). 

Part (b). Suppose that X satisfies (28). Then the sets Bi defined in (40) 
satisfy v{Bi) > for % — 1, . . . , d. Consequently, the random variables \X®\ 
are regularly varying with tail index a. Denoting 

«;« := (i/(fl i ))~ 1/a , (51) 
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> t/w®} 


P{|X«: 


>t} 



one obtains 

P{\w®X®\>t} 

£JS> pqxw^t} = t^y p{|x«|>t} p{||x|| 1 >t} 

= 1 

for i — 1, . . . , d. Hence, for any i,j e {1, . . . , d}, 

P{\w®X®\ > t\ 

lim ^ TTT 7" = 1. 

t^oo P{| W 0)X0)| > i} 

To prove the inverse implication, suppose that Z := wX has balanced 
tails for some w G (0, oo) d . Then the the exponent measure v of X satisfies 

v(Bj) _ P{\xM\>t} 
u{B 1 ) t™ P{|X«| >t} 

P{|Z«| > w^t} 



= lim 



t^ooP{\ZW\ > wWt} 

fp{\Z^\>w^t) 
t 1 ™ I P{|Z«|>t} 







>t} P{\ 




>t} 




>w( 1 H} P{\ZW 


\>t} 



W 



(1) 



e(0,oo), ie {!,..., d} 



Since multivariate regular variation of X implies v(Bj) > for at least one 
index j 6 {1, . . . , d}, this yields i^(-Bj) > for all i. 

Part (c). This is an immediate consequence of part (a) and the invariance 
of canonical spectral measures under componentwise rescaling. □ 

Representation (49) suggests that ordering of the normalized extreme risk 
indices 7^/(7^ + 7_ ei ) in the balanced tails setting can be considered as an 
integral order relation for canonical spectral measures with respect to the 
function class 

G a ■■ {<k.„-i( ^j- (52) 
This justifies the following definition. 

Definition 3.4. Let \P* and $* be canonical spectral measures on Sf and let 
a > 0. Then the order relation fy* ^g a $* is defined by 

VgeG a y*g<$*g. (53) 
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Remark 3.5. (a) For a = 1 and spectral measures on S d the extreme risk in- 
dex 7£pT) is linear in £ (cf. Mainik and Riischendorf, 2010, Lemma 3.2). 
Consequently, ^g a is indifferent in this case, i.e., any \&* and $* on £>(£ d ) 
satisfy 

^* =<0,i $* and $* r< 0)1 (54) 

(b) The order relation is mixing invariant in the sense that uniform 
ordering of two parametric families {f J : i5 G 0} and {$J : $ G 6}, 

implies 



./e 

for any probability measure /x on 0. 

The following theorem states that ^ ap i is in a certain sense equivalent to 
the ordering of canonical spectral measures and allows to reduce the verifica- 
tion of ^ ap i to the verification of ^g a . Some exemplary applications are given 
in Section 5. Furthermore, given explicit representations of spectral measures 
or their canonical versions, this result allows to verify ^ ap i numerically, which 
is very useful in practice. 

Theorem 3.6. Let X andY be multivariate regularly varying random vectors 
on M. d with tail index a G (0, oo) and canonical spectral measures ty* x and 
\&y. Further, suppose that X andY satisfy the balanced tails condition (48). 

(a) If ^ apl |yW|, then ** x < Qa implies X ^ apl Y. 

(b) If\X<V\ ^ apl and ^ apl then^* x <g a is equivalent 
toX^Y. 

Proof, (a) Since X has balanced tails, Proposition 3.3(a) yields 

p {ex > t} f p{ex>t} p{\\x\\ 1 >t} 

ESoP{\XW\>t} «™ ^P{||X||i >*} 'p{\xw\ >t} 

jdx) 

7ei (X) + 7 - ei (X) 



Analogously one obtains 



p\ey > t\ 

}^ P{\Y^\>t} = ^ 9 ^ 
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Moreover, ^f* x <g a ^ Y implies 



< 1. (55) 



Consequently, 



= lim sup 

< 1 







\>t} 


p{|yw| 


>t} 







>;} 


p{|y(*)| 


>t} 



lim sup (56) 



due to (55) and ^ api |y«|. 

(b) By part (a), it suffices to show that X ;< ap i Y implies ^>* x -<g a By 
assumption l-X^ 1 )] and \Y^\ have asymptotically equivalent tails, 

P{\XW\ >t} 



lim 



Thus (56) yields 



t^oc P{|y(i)| > t} 
- hmsup 



and X ^ apl y implies ^ ^ Q V* Y . □ 

The following result answers the question for dependence structures cor- 
responding to the best and the worst possible diversification effects for mul- 
tivariate regularly varying random vectors in According to Theorem 3.6, 
it suffices to find the upper and the lower elements with respect to <g a in the 
set of all canonical spectral measures on S d . It turns out that for a > 1 the 
best diversification effects are obtained in case of asymptotic independence, 
i.e., the ^g Q -maximal element is given by 

d 

i=i 

whereas the worst diversification effects are obtained in case of the asymptotic 
comonotonicity, represented by 

n-=d-8 [l/d _ l/d) . (58) 

For a < 1 the situation is inverse. 
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Theorem 3.7. Let ^* be an arbitrary canonical spectral measure on S d and 
let y&g and ^ be defined according to (57) and (58). TTien 

(a) % ±g a ** ±g a nfora>l. 

(b) m\ ±g a V* ±g a % for a G (0,1]. 

Proof. Let X be multivariate regularly varying on with canonical spectral 
measure Without loss of generality we can assume that X satisfies the 
balanced tails condition (48). Then, according to (49), we have 

= m 

Furthermore, we have ^*g eua = 1 for % — 1, . . . , d. Recall that the mapping 
£ i— > 7,t is convex for a > 1 (cf. Mainik and Riischendorf, 2010, Lemma 3.2). 
Due to (59) this behaviour is inherited by the mapping £ h-> ^g^a- Thus 
for a > 1 we have ^*g^, a < 1 = **^,a for all £ G S d , which exactly means 
^* *I for a > 1. 

To complete the proof of part (a), note that the normalization of canonical 
spectral measures yields 



i=l ■ /sd i=l 



(60) 



Comparing the integrand on the right side of (60) with the function g^ a ( s ) 
(£ T s 1/a ) a , we see that 



d 

a 



i=l i=l 

with 

f« . ( s W) 1/a 

Thus it suffices to demonstrate that Ylt=i z ? — x ) which follows from Zi G 
[0, 1], zf < Zi for a > 1, and E^ =1 ^ = 1. 

The inverse result for a G (0, 1] stated in (b) follows from the concavity 
of the mapping £ \-t ^*g^ a and the inequality zf > z\. □ 

Due to Theorem 3.6, an analogue of the foregoing result for ^ ap i is 
straightforward. 
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Corollary 3.8. Let X be multivariate regularly varying in M.+ with tail in- 
dex a G (0, oo) and identically distributed margins ^ F , i = 1, . . . , d. 
Further, let Y be a random vector with independent margins Y® ~ F, and 
let Z be a random vector with totally dependent margins 

Z {i) = Z (i) p_ a _ s _ 

and Z^ ~ F . Then 

(a) Y ^ ap i X ^ apl Z for a > 1 

(b) Z^X^Y for a e (0,1]. 

Remark 3.9. The strict assumptions of Corollary 3.8 are chosen for clear- 
ness and simplicity. The independence of and the total dependence of 
are needed only in the tail region, i.e., it suffices for Y and Z to be 
multivariate regularly varying with canonical spectral measures \1/q and ^l, 
respectively. Furthermore, the assumption of identically distributed margins 
can be replaced by equivalent tails: 

, P{F« >t\ , r ;z > / ; 
1 = lim — ; — = lim — £ — y-, i — 1, . . . , a. 



p {r« > t\ p {z^ > t\ 

= lim — ; — ttt 7- = lim — ^ — ttt 

t-^oo P {X« > *} t^oo P {X« > *} ' 



Finally, the non-negativity of X^\ Y^ l \ and is needed only in the asymp- 
totic sense. The ordering results remain true if the spectral measure of X is 
restricted to the unit simplex S d . 

Combining Theorem 3.6 with Theorem 2.7, one obtains an ordering re- 
sult for the canonical spectral measures of multivariate regularly varying el- 
liptical distributions. The notation = ^/*(a,C) is justified by the fact 
that spectral measures of elliptical distributions depend only on the tail 
index a and the generalized covariance matrix C. An explicit representa- 
tion of spectral densities for bivariate elliptical distributions was obtained by 
Hult and Lindskog (2002). Alternative representations that are valid for all 
dimensions d > 2 are given in Mainik (2010), Lemma 2.8. 

Proposition 3.10. Let C and D be d-dimensional covariance matrices sat- 
isfying 

C M = A,i > 0, i = l,...,d, (61) 

and 

V£e£ d £ T C£<£ T £>£. (62) 

Then 

V«>0 V*(a,C)lg a S!*(a,D). 
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Proof. Fix a G (0, oo) and consider random vectors 

X = RAU, Y = RBU, 

where A and B are square roots of the matrices C and D in (62), i.e., 

C = AA T , D = BB T , 

and R is an arbitrary regularly varying non-negative random variable with 
tail index a. 

As a consequence of Theorem 2.7 one obtains X ;< ap i Y . Furthermore, 
invariance of ^ ap i under componentwise rescaling yields wX ;< ap i wY for 
w = (w {1) ,...,w id) ) with 

w 0) : = C M - 1/2 = A,i" 1/2 , i = l,...,d. 

Moreover, as a particular consequence of arguments underlying (22), one 
obtains 

w ®x® ±w {j) Y {j \ i,j e {l,...,d}. 

Hence the random vectors wX and wY satisfy the balanced tails condition 
(48), whereas their components are mutually ordered with respect to z< ap i. 
Finally, Theorem 3.6(b) and invariance of canonical spectral measures under 
componentwise rescalings yield 

C) = Kx ±Q a Ky = D). □ 



The subsequent result extends Theorem 3.6 to random vectors that do 
not have balanced tails. 

Theorem 3.11. Let X and Y be multivariate regularly varying random vec- 
tors on M. d with tail index a G (0, oo) and canonical spectral measures 
and \l/y. Further, assume that ^ ap i \Y^\ with 



P{l 




\>t} 


p{|y(*)| 


>t} 



X i ■= J im . D rv m >t J e (°> 1] ( 63 ) 
for i = 1, . . . , d and that the vector v = (v^\ . . . , v^) defined by 

u« := \; 1/a (64) 

satisfies 

X^vX or v^Y^iY. (65) 
Then V* x ^g a ^ implies X ^ api Y. 
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Proof. According to Proposition 3.3(b), there exists w G IR+ such that wY 
satisfies the balanced tails condition (48). Furthermore, the tails of the ran- 
dom vector 

vwX := (v^w^X™, . . .^W")) 
with v defined in (63) are also balanced. Indeed, it is easy to see that 

p/UWyWI > A P{\v®w®X®\ > t\ , ,. UQ 

lim — \} r 1 — = lim — ^ F . ... 1 — r-^- = u> w 

for i — 1, . . . , d. Analogously one obtains 

P{pW)|>t} ^ 
*™ P{|XW|>t} 1 J 

and, as a result, 
lim 



p{l 




>*} 


p{| w »y«) 


><} 



lim 
lim 



P{l 




>;} 


p{|y(0| 


>0 



P{| 




>;} 


p{|a«: 









\>t} 


p{|y«| 


>t} 



P{|X«|>t} 

= Ai '/™ p{|yw|>0 
= l 

for i = 1, . . . , d. Hence the balanced tails condition for wY implies that the 
tails of vwX are also balanced. 

Furthermore, invariance of canonical spectral measures under component- 
wise rescaling yields 

K W X = ±G a = Ky- 
Thus, applying Theorem 3.6(a), one obtains 

vwX ^ ap i wY. (66) 
Since = \ > for % = 1, . . . , d, condition (66) is equivalent to 

wX ^ apl v-'wY. (67) 
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Moreover, assumption (65) implies 

wX ^ ap i vwX or v ~ l wY < &v \ wY. (68) 
Combining this ordering statement with (66) and (67), one obtains 

WX ^ ap i WY. 

Finally, invariance of ^ ap i with respect to componentwise rescaling yields 
X^Y. ' " □ 

In the special case of random vectors in Theorem 3.11 can be simplified 
to the following result. 

Corollary 3.12. Let X and Y be multivariate regularly varying random 
vectors on IR+ with tail index a G (0, oo) and canonical spectral measures 
and Further, suppose that 

P ( \X { - €) \ > t\ 

K :-hmsup 1 { G(0,1], i = l,...,d. (69) 

Then < Qa m* Y implies X ^ ap i Y. 

Proof. Assumption (69) yields that the rescaling vector v defined in (64) is 
an element of [1, oo) d . Thus v — (1, . . . , 1) G and, since X takes values in 
R+, we have 

X^ pl X + (v-(l,...,l))X = vX. 

Similar arguments yield v~ x Y Y. Hence condition (65) of Theorem 3.11 
is satisfied. □ 

The final result of this section is due to the indifference of -<g a for a — 1 
mentioned in Remark 3.5(a). This special property of spectral measures on 
S d allows to reduce ^ ap i to the ordering of components. It should be noted 
that this result cannot be extended to the general case of spectral measures 
on Sf. 

Lemma 3.13. Let X and Y be multivariate regularly varying on with 
tail index a — 1. Further, suppose that Y satisfies the non- degeneracy con- 
dition (28) and that X« ^ api for i = 1, . . . , d. Then X ^ apl Y. 

Proof. According to Proposition 3.3(b), there exists w G (0, oo) d such that 
wY satisfies the balanced tails condition (48). Furthermore, due to the in- 
variance of ^ ap i under componentwise rescaling, X ^ ap i Y is equivalent to 
wX ^ ap i wY. 
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Thus it can be assumed without loss of generality that Y has balanced 
tails. This yields 

P {X« >t\ P {X« > t\ 

Ai := U f-L up p {r« >t}= ^ p { yd) > o ' < = 1 '-' d - 

Hence the assumption ;< ap i for i = 1, . . . , d implies Aj G [0, 1] for all 
i. Moreover, the balanced tails condition for Y yields 

7ei on = ... = 7 e d (n (to) 

Now consider the random vector X and denote 

j := argmax7 6i (X). 

ie{l,...,d} 

Recall that 7e i (- 5 — vx{{% G 1}) with ux denoting the exponent 

measure of X and that vx is non-zero. This yields 7e .(X) > even if X 
does not satisfy the non- degeneracy condition (28). Moreover, for a — 1, the 
mapping £ i— >■ 7g(X) is linear. This implies 

d 

^( x ) = • ^ e G s " ( 71 ) 

and (70) yields 

7^ = ££ W -7 ei (n = 7 ei (n (72) 
i=i 

Hence 

limsup^^i 

TT P p {£ T r > t} 

fp{CX>t} P{X®>t} P{Y^>t} 
~ 1 ™ S " P ^P{X0-)>t} ' P{F(D>t} ' P{CY >t} 

- 7ej (X) J 7f (F) 
< 1 

due to Aj < 1, (71), and (72). □ 
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4 Relations to convex and supermodular or- 
ders 



As mentioned in Remark 2.4(b), dependence orders ^ sm , ^dcx and convexity 
orders ^ cx , z^i CX! Implex do not imply ;< ap i in general. However, it turns out 
that the relationship between ^ apl and the ordering of canonical spectral 
measures by ^g a allows to draw conclusions of this type in the special case 
of multivariate regularly varying models. The core result of this section is 
stated in Theorem 4.1. It entails a collection of sufficient criteria for ^ ap i in 
terms of convex and supermodular order relations, with particular interest 
paid to the inversion of diversification effects for a < 1. An application to 
copula based models is given in Proposition 4.4. 

This approach was applied by Embrechts et al. (2009b) to the ordering of 
risks for the portfolio vector £ = (1, . . . , 1) and for a specific family of mul- 
tivariate regularly varying models with identically distributed, non-negative 
margins A"W ( c f. Example 5.2 in Section 5). 

The next theorem is the core element of this section. It generalizes the 
arguments of Embrechts et al. (2009b) to multivariate regularly varying ran- 
dom vectors in M. d with balanced tails and tail index a^l. The case a = 1 
is not included for two reasons. First, this case is partly trivial due to the 
indifference of -<g a for spectral measures on E d (cf. Remark 3.5(a)). Second, 
Karamata's theorem used in the proof of the integrable case a > 1 does not 
yield the desired result for random variables with tail index a = 1. 

Theorem 4.1. Let X and Y be multivariate regularly varying on M. d with 
identical tail index q^I. Further, assume that X and Y satisfy the balanced 
tails condition (48). 

(a) For a > 1 let 

p{\xm\>t\ , x 

l 7J^p{\YV\>t } = 1 (73) 

and let there exist u > such that with h u (t) := (t — u) + 

Vu > u V£ e H d Eh u (£ T X) < Eh u (£ T F) . (74) 
Then ** x <g a ^* Y . 

(b) For a < 1 suppose that \X^\ and \Y^\ are equivalent with respect to 

\X^\ ^ apl |yW| and |yW| ^ apl |X«|, (75) 
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and let there exist u > such that with f u (t) := — (t A u), 

Vu>u vtez d E/ u ((e T x) + ) <E/ u ((e T r) + ). (76) 

Then y* Y ± Qa V* x . 

The proof will be given after some conclusions and remarks. In particular, 
it should be noted that the relation between -<g a and ;< ap i established in 
Theorem 3.6 immediately yields the following result. 

Corollary 4.2. (a) If random vectors X and Y satisfy conditions of Theo- 
rem 4.1(a), then X ^ ap i Y; 

(b) If X and Y satisfy conditions of Theorem \.l(b), then Y ^ ap i X . 

It should also be noted that conditions (74) and (76) are asymptotic forms 
of the increasing convex ordering ^ X ^ icx £ T F and the decreasing convex 
ordering £ T X z<dc CX £ T ^, respectively. The consequences can be outlined as 
follows. 

Remark 4.3. (a) The following criteria are sufficient for (74) and (76) to 
hold: 

(i) (£ T X) + ^ cx (CY) + for all f G £ d , 

(ii) X and Y are restricted to and X ■< Y with ■< denoting either 

Implex; I^lcx? — cxj ZMcxj Or I^ sm . 

(b) Additionally, condition (74) follows from X ^ Y with ^ denoting either 

Implex; I^lcx; — Cxi I^dcx; Or Z^ sm . 

Finally, a comment should be made upon convex ordering of non-integrable 
random variables and diversification for a < 1. The so-called phase change 
at a — 1, i.e., the inversion of diversification effects taking place when the 
tail index a crosses this critical value, demonstrates that the implications 
of convex ordering are essentially different for integrable and non-integrable 
random variables. Indeed, it is easy to see that if a random variable Z on R 
satisfies E[Z + ] = E[Z_] = 00, then the only integrable convex functions of Z 
are the constant ones. Moreover, if Z is restricted to M + and EZ = 00, then 
any integrable convex function of Z is necessarily non- increasing. 
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Proof of Theorem ^.i.(a) Consider the expectations in (74). It is easy to see 
that for u > 

-Eh u (Cx) = - [ P{CX>t}dt 
= I ¥{i T X>tu}dt 

7(1,00) 

and, as a consequence, 

u^Ehu (£ T X) P {£ T X > u} f P {£ T X > tu] 



p{|xw|>«} p{|x(D|>«}y (1>00) p{e T x> M } dl 

Moreover, Proposition 3.3(a) implies 

il n i P{|.Y<.»|>„. } = 7tl (A-) + 7- ei (.Y) = (77) 

and Karamata's theorem (cf. de Haan and Ferreira, 2006, Theorem B.1.5) 
yields 

r p\Cx>tu\ r i 
iim / — f dt = / r a dt = — -. 

™7(l,oo) P{£ T X>«} 7(1,00) 

As a result one obtains 

.. nr^EK (Cx) 1 

lm = Wyflt„ 

«^oo P{|X( 1 )| > «} a-1 



and, analogously, 



lim 



u- l Eh u (£ T y) 



«->oo P{|y(i)| >u} a -l^* Y9( ' a ' 
Hence (74) and (73) yield 

u-'ek (ex) 



lim sup 

u— >oo 



u 



- i Eh u (ex) p{|yw|> u } p{\x^\> u y 



P{\XW\>u} u- l Eh u (CY) P{\YW\>u} 



for all £ G S d , which exactly means ty* x <g a *y. 
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(b) Note that (75) implies 

p{\xw\>t\ , x 

^ pTO 1 (78) 

and that (76) yields 

Vu > «<, Vu > E/ u+t , (e T X) - E/ tt+w {Cy) < 0. (79) 
Furthermore, it is easy to see that any random variable Z in R + satisfies 

E[Z Au] = [ (tAu)dP z (t) 

</(0,oo) 

= / / l{s < t} ■ l{s < u}dsdP z (t) 

i(0,oo) </(0,oo) 

= / l{s<u} l{s < t}dP z {t)ds 

J(0,oo) J(Q.oo) 



'(0,oo) ^(0,oo) 

f 

'(0,u) 



[ P{Z>s}ds. 

J(0,u) 



This implies 



Ef u+v (Z) = Ef u (Z) - ! P{Z>t} dt. 

J (u,u+v) 



Consequently, (79) yields 

\/u > u Vv > Ef u (CX) - Ef u (£ T Y) < I(u, v) (80) 



where 



I(u, v) := I (P {CX >t}-P {^Y > t}) 

J (u,u+v) 

= [ <P(t)-P{\X^\>t}dt 

J (u,u+v) 



with 



P \£ T X >t\ -P \£ T Y >t\ 



P{\X0-)\>t} 

Moreover, (78), (77), and an analogue of (77) for Y yield 

P{£ T X>t} P{CY>t} P{\YW\>t} 



0(t) 



P{\XW\>t} P{\YW\>t} P{\XW\>t} 
-)• V* x gt, a - y* Y gt, a , t OO. (81) 
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Now suppose that \l/y -<g a is not satisfied, i.e., there exists £ G S d 
such that \&y<7£, a > ^fl^a- Then (81) yields <f>(t) < —e for some e > and 
sufficiently large i. This implies 

I(u,v)<-e I P > t) dt (82) 

J (u,u+v) 

for sufficiently large u and all t> > 0. Moreover, regular variation of 

with tail index a < 1 implies E = oo. Consequently, the integral on 

the right side of (82) tends to infinity for v — > oo: 

Vw > lim / P { |X (1) I > t} dt = oo. 

v ^°° J(u,u+v) 

Hence, choosing u and v sufficiently large, one can achieve I(u,v) < c for 
any cel. In particular, u and v can be chosen such that 

I(u,v)<Ef u (ex)-Ef u (eY), 

which contradicts (80). Thus ^yg^ a > ^* x 9(,,a cannot be true and therefore 
it necessarily holds that \I/y -<g a ^* x . □ 

Now let us return to the ordering criterion in terms of the supermodular 
order ;< sm stated in Remark 4.3. The invariance of ^ sm under non-decreasing 
component transformations allows to transfer these criteria to copula models. 
Furthermore, since we are interested in the ordering of the asymptotic de- 
pendence structures represented by the canonical spectral measures, \&* and 

we can take any copulas that yield and ^ as asymptotic dependence 
structures. 

A natural choice is given by the extreme value copulas, defined as the 
copulas of simple max-stable distributions corresponding to i.e., the dis- 
tributions 

G*{x) :— exp (—u* (-[oo,x] c )) , x eR d + (83) 

where v* is the canonical exponent associated with \1>* via (41). For further 
details on max-stable and simple max-stable distributions we refer to Resnick 
(1987). Since extreme value copulas and canonical spectral measures can be 
considered as alternative parametrizations of the same asymptotic depen- 
dence structures, we obtain the following result. 

Proposition 4.4. Let and ^ be canonical spectral measures on H d . Fur- 
ther, fori = 1,2, let C% denote the copula of the simple max-stable distribution 
G* induced by \&* according to (83) and (41). Then C\ ^ sm C2 implies 
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(a) n ^ a f5/orae(l,oo); 



(b) < Ga *J for a e (0,1). 

Proof. Let z/* denote the canonical exponent measures corresponding to \&* 
and G*. It is easy to see that the transformed measures 

v a ,i-.= v* oT~ l , i = l,2, 
with a > and the transformation T defined as 

T-. x ^ ((x«) 1/Q ,...,(^) Va ), zeRj, 
exhibit the scaling property with index —a: 

i/ Q)i (L4) = rv aii (A), A e B(R d + \ {0}). 
Hence the transformed distributions 

G a ,i(x) :=G*oT- 1 (x) = exp(-u^([0,x] c )) (84) 

are max-stable with exponent measures v a ^. 

It is well known that max-stable distributions with identical heavy-tailed 
margins are multivariate regularly varying (cf. Resnick, 1987). Moreover, 
the limit measure v in the multivariate regular variation condition can be 
chosen equal to the exponential measure associated with the property of 
max-stability. Consequently, the probability distributions G a ^ for i = 1,2 
and a > are multivariate regularly varying with tail index a and canonical 
spectral measures 

Furthermore, it is easy to see that X ~ G ail and Y ~ G a> 2 have identical 
margins: 

X«4yw ije{l,...,d}. 

Moreover, due to the invariance of ^ sm under non-decreasing marginal trans- 
formations, C\ ^ sm C2 implies 

G a ,l ^sm G a> 2 

for all a > 0. Thus an application of the ordering criteria from Remark 4.3 
to X ~ G at i and Y ~ G a> 2 completes the proof. □ 
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5 Examples 

This section concludes the paper by a series of examples with parametric 
models illustrating the results from the foregoing sections. Examples 5.1 and 
5.2 demonstrate application of Proposition 4.4 to copula based models and 
the phenomenon of phase change for random vectors in The fact that the 
phase change does not necessarily occur in the general case is demonstrated 
by multivariate Student-t distributions in Example 5.3. 

Example 5.1. Recall the family of Gumbel copulas given by 



Gumbel copulas are extreme value copulas, i.e., they are copulas of simple 
max-stable distributions. According to Wei and Hu (2002), Gumbel copulas 
with dependence parameter e [1, oo) are ordered by ^ sm : 



Consequently, Proposition 4.4 applies to the family of canonical spectral 
measures corresponding to the Gumbel copulas C#. Thus 1 < 0i < 02 < 
00 implies -<g a \f^ 2 for a > 1 and there is a phase change when a crosses 
the value 1, i.e., for a e (0, 1) there holds \I>Jj 2 ^g a 

Applying Theorem 3.6, one obtains ordering with respect to ^ ap i for ran- 
dom vectors X and Y on IR^_ that are multivariate regularly varying with 
canonical spectral measures of Gumbel type and have balanced tails ordered 
by ^api- In particular, this is the case if X and Y have identical regularly 
varying marginal distributions and Archimedean copulas that satisfy appro- 
priate regularity conditions (cf. Genest and Rivest, 1989; Barbe et al., 2006). 

Moreover, it is also worth a remark that multivariate regularly varying 
random vectors with Archimedean copulas can only induce extreme value 
copulas of Gumbel type (cf. Genest and Rivest, 1989). 

Figure 1 illustrates the resulting diversification effects in the bivariate 
case, including indifference to portfolio diversification for a = 1 and the 
phase change occurring when a crosses this critical value. The graphics show 
the function £^ t-> ^g^ a for selected values of and a. Due to X e M+, 
representation ^g^ a = 7?/(7 ei + 7-eJ simplifies to ^%g^ a = l^/la and 
therefore 




(85) 



V0i, 02 e [l, 00) 0i < 2 ^ CV 



(86) 



^<?ei,a = %9, 



= 1. 
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Extreme risk index (normalized) 



Extreme risk index (normalized) 





(b) Varying a for -Q = 2 
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(c) Varying $ for a > 1 



(d) Varying $ for a < 1 



Figure 1: Bivariate Gumbel copulas: Diversification effects represented by 
functions £W >->■ ^fi^.a for selected values of $ and a. 
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As already mentioned above, Theorem 4.1 generalizes some arguments 
from Embrechts et al. (2009b). The next example concerns Galambos copu- 
las as addressed in that original publication. 

Example 5.2. Another family of extreme value copulas that are ordered by 
^ sm is the family of d-dimensional Galambos copulas with parameter i? G 
(0,oo): 

C7,(u) := exp ( £ (-1)' 1 ' ( £ (- logii«) ]. (87) 

\/c{i,...,d} \«e/ J J 

According to Wei and Hu (2002), -&\ < i? 2 implies C.g 1 ^ sm C$ 2 . Thus 
Proposition 4.4 yields ordering of the corresponding canonical spectral mea- 
sures ^ with respect to ^g a . Similarly to the case of Gumbel copulas, 
A < $ 2 implies ±g a % 2 for a > 1 and ^ ^ for a G (0, 1). 

Finally, it should be noted that Galambos copulas correspond to the 
canonical exponent measures of random vectors X in with identically 
distributed regularly varying margins X® and dependence structure of — X 
given by an Archimedean copula with a regularly varying generator 0(1 — 
1/t). Models of this type were discussed in recent studies of aggregation 
effects for extreme risks (cf. Alink et al., 2004, 2005; Neslehova et al., 2006; 
Barbe et al., 2006; Embrechts et al., 2009a,b). 

The final example illustrates results established in Proposition 3.10 and 
Theorem 2.7. In particular, it shows that elliptical distributions do not ex- 
hibit a phase change at a = 1. 

Example 5.3. Recall multivariate Student-t distributions and consider the 
case with equal degrees of freedom, i.e., 

X = fi X + RA X U, Y = fi Y + RA Y U, (88) 

where R = \Z\ for a Student-t distributed random variable Z with degrees of 
freedom equal to a G (0, oo). Further, let the generalized covariance matrices 
Cx = C(px) and Cy = C(py) be defined as 

C( P ) ■= ( I J ) (89) 

and assume that px < Py- 

As already mentioned in Remark 2.8(a), Cx and Cy satisfy condition (62) 
and Proposition 3.10 yields X ^ ap i Y. Moreover, Proposition 3.10 implies a 
uniform ordering of diversification effects in the sense that 

n = K, Px -<g a K, PY = 
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for all a E (0, oo). 

Figure 2 shows functions £^ >->■ ^p^a f° r selected parameter values p 
and a that illustrate the ordering of asymptotic portfolio losses by p and the 
missing phase change at a — 1. The indifference to portfolio diversification 
for a = 1 is also absent. Moreover, symmetry of elliptical distributions 
implies 7_ ei = j ei and, as a result, 

Thus the standardization of the plots in Figure 2 is different from that in 
Figure 1. 

Remark 5.4. All examples the authors are aware of suggest that the diversi- 
fication coefficient ^*g^ a is decreasing in a. This means that risk diversifi- 
cation is stronger for lighter component tails than for heavier ones. 

However, it should be noted that the influence of the tail index a on risk 
aggregation is different from that. The asymptotic risk aggregation coefficient 

P{lW + ... + lM >t\ 
:= lim — — -—, r - 

HA t^oo P {XW > t} 

introduced by Wuthrich (2003) is known to be increasing in a when the loss 
components are non-negative (cf. Barbe et al., 2006). It is easy to see 
that the restriction to non-negative implies 

,. P{\\Xh>t} 1 

qd -}^ p {x w>t} - 

Moreover, denoting the uniformly diversified portfolio by rj, 

V^d-^l,...,!), 

one obtains 

,. P WX > d-H) lv 

q d = hm — j- - r — ttt -r- 1 - = a — -. 

H t^oc P {XW > t} 7ei 

Thus qa is a product of the factor d a , which is increasing in a, and the ratio 
7r?/7ei, which is closely related to to the diversification coefficient ^*g^ a . 
In particular, given equal marginal weights, i.e., 

7ei — ••• — 7e d ? 

Proposition 3.3(a) yields 

7ei 
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Extreme risk index (normalized) 



Extreme risk index (normalized) 
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Figure 2: Bivariate elliptical distributions with generalized covariance ma- 
trices defined in (89): Diversification effects represented by functions (f*- 1 - 1 i—>- 
^*a,p9(,,a f° r selected values of p and a. 
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As already mentioned above, the coefficients ^*g^, a with (GE d are decreas- 
ing in all examples considered here. This means that the aggregation and the 
diversification of risks are influenced by the tail index a in different, maybe 
even always contrary ways. 

The question for the generality of this contrary influence is currently open. 
One can easily prove that the extreme risk index 7^ = ^f^ a is decreasing 
in a for £ e S d . However, this result cannot be extended to ^*g^ a directly 
since ty*g^ a is related to ^f^, a by the normalizations (49) and (50). The 
question whether ty*g^ a with arbitrary £ G E d or at least ^*g v ^ a is generally 
decreasing in a is an interesting subject for further research. 
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